Controllability of NLS in the vicinity of solitary 

wave solutions 

H. Lange (Cologne) 
H. Teismann (Acadia) 

February 2, 2008 

Abstract 

Local exact controllability of the ID NLS (subject to zero boundary con- 
ditions) with distributed control is shown to hold in a //^-neighbourhood of 
the nonlinear ground state. The Hilbert Uniqueness Method (HUM), due to 
J.-L. Lions, is applied to the linear control problem that arises by linearization 
around the ground state. The application of HUM crucially depends on the 
spectral properties of the linearized NLS operator which are given in detail. 

1 Introduction 

The control properties of many PDEs arising in physics and engineering have been 
studied extensively. Those investigations include exact and/or optimal controlla- 
bility for the (linear and nonlinear) heat, wave, beam and plate equations as well 
as equations of elasticity and Navier Stokes equations, to mention just some of the 
most prominent examples. For Schrodinger equations, however, the control theory 
is markedly less developed. For a brief survey on control results for (linear and non- 
linear) Schrodinger equations, see e.g. J7j. The purpose of this paper is to establish 
local exact controllability for the ID cubic nonlinear Schrodinger equation (subject 
to zero boundary conditions) with distributed control. Specifically, we will consider 
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the following control problem. 



(la) iy t 

(lb) y(t,0) 

(lc) y(0,x) 

(Id) y(T,x) 



-y xx + f{\y\ 2 )y + g al {x)<t^), e (o, l), t e (o,T)) 
y(U) = (t€[0,T]) 
2/ (x) (x 6(0,1)) 
yi(x) (x 6(0,1)) 



where g n ,( x ) denotes the indicator function for some fixed, possibly small, open 
subinterval Q' C Q := (0, 1) which represents the spatial region in which the control 
is applied. Given a fixed control "horizon" T > and initial and target states 2/0 
and j/i, the objective is to construct a control function u = u(t,x) that will steer 
the state y from y to 2/i> i-e. the unique solution y = y(t,x) of (|Ta j) -(fTc |) is to 
satisfy (fid)). The control problem (|laj ) -(|ld j) was posed in [T71 [TB]. A small-data 
controllability result for periodic boundary conditions is contained in 
In the remainder of this paper, we will concentrate on the special case of the focusing 
cubic nonlinear Schrodinger equation (NLS), i.e. 



More general nonlinearities could be treated, but this is not our main interest here. 
So we restrict ourselves to the prototypical cubic nonlinearity. Our main result states 
that the control problem (jlaj) - (jld|) is soluble locally within an if ^neighbourhood 
of the (ground-state) solitary-wave (or soliton) solution ip^t, x) = e* M '0 M (x), if the 
control time T > is sufficiently large. Here <p^(x) denotes the (nonlinear) ground 
state 1 , i.e. the (real and) positive solution of the boundary value problem 



which is known to exist and to be unique (see Section lA~Tf . Our main results reads. 

Theorem 1 Let Q' C (0,1), fi > be given and let M denote the ground state. 
Then there exist T > and 6 > such that, for any y , y\ 6 Hq(0,1) := 
{v 6 H\0, 1) | v (0) = v(l) = 0} satisfying 



= —s. 



(2a) 
(2b) 



- 4>" + fx(f> - <j) 3 = o (x 6 n) 

0(0) = 0(1) = o. 



Vo - <t>v\\m < $ and 1 1 2/1 - ^\\ H i < 5, 



Often the time-dependent solution ip^it^x) is referred to as the ground state. 
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the control problem Ala)) - I7dj) has a solution u = u(t,x), i.e. there exists a control 
function u G C([0, T}; Hq(0, 1)) such that the unique solution y G C([0,T]; Hq (0, 1)) 
of < f7a)) -< f7c)) satisfies I7dj) . 

The theorem will be proved by applying the implicit function theorem (IFT) to the 
nonlinear map ^ : [H%(0, l)] 2 x C([0, T];Hq(0, 1)) -> #■}((), 1), defined by 

(3) *(2/o,2h,«) := y(T;y ,u) - y x , 

where t i— > y(t;yo,u) G C([0, T]; ifg(0, 1)) denotes the unique solution of (fTa j) -([Tc |) . 
To be able to apply the IFT, it will be verified that the linearization <9 U \I/(0 M , e ttlT 4>n, 0) : 
C([0, T]; ^(0, 1)) -> #*((), 1) of at the point e^ T M , 0) exists as a bounded 
map and possesses a bounded inverse. This amounts to showing that the linear PDE 
that arises by linearizing NLS around the stationary solution tp^ is exactly control- 
lable; this, in turn, is done by employing the Hilbert Uniqueness Method (HUM) due 
to J.-L. Lions 0. The main difficulty in the application of HUM stems from the lack 
of selfadjointness of the linearized operator, which makes determining its spectral 
properties more intricate. The analysis reveals that most spectral properties known 
to hold for the linearized NLS operator in the whole-space case |1(J[ H], carry over 
to the zero-boundary case considered in this paper (see Section 13. 1.1JI 2 . This fact is 
of independent interest, but it is, to the best of our knowledge, not available in the 
literature. 

This paper is organized as follows. In Section 2 we formulate the linear control 
problem that arises by linearization around the ground state (system (J9aj) - (j9d|) ) and 
state its solvability (Th. |2J). We also show how Theorem ^ is derived from this 
controllability result ( Section 12. 2j) . Section 3 contains the proof of Th. The proof 
is based on HUM and, thus, hinges on the "observability estimates" that are proved 
in Sections 13.11 (L 2 -estimate) and 13.21 (if ^estimate). The all- important spectral 
properties needed in these proofs are listed in Section 13.1.11 and proved in the Ap- 
pendix (|B.3|) . The Appendix also contains information on explicit solution formulas 
for (J2"aj) . (f2b"|) f Section |A.1|) . a variational description of the ground state (jA.2|) . and 
asymptotic formulas for the eigenvalues and eigenfunctions (jB.2|) . A number of open 
problems are listed in Section 0] 

2 Obviously, this statement is not meant to be applied to those properties which crucially depend 
on the fact that the spatial domain f2 = (0, 1) is bounded, such as the absence of a continuous part 
of the spectrum. 
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2 Linearization and proof of Theorem U 



Let T > be such that property ^151 in I3.1.1l is satisfied. (Such a T exists according 
to ^4H1 in IB~3l ) The parameters T > 0, /i > and the interval Q' C (0, 1) will be 
kept fixed in all what follows. We also assume that g := g n , is a smooth function 
satisfying supp(g) C Q! and < g(x) < l. 3 

2.1 Linearization of \£ 

We write \I/ as 

*(yo,yi,«) = $(yo,«) - j/i, 

where $ : ^(0, 1) x C([0, T]; H^O, 1)) -> ^(0, 1) is the map $(y , u) := y , 
To see that the map $ (and hence ^) is well-defined, we need to know that the initial 
value problem (|Ta | -(fIc jl has a unique solution y 6 C([0,T]; i^o(0, 1)) for any choice 
of data y e #o(0, 1) and «6 C([0,T]; Hq(0,1)). This is known for the homogeneous 
NLS (i.e. it = in (jla|) ) in ID; cf. |3J Corollary 3.5.2.]. It is fairly easy to convince 
oneself that NLS with an additional inhomogeneity u G C([0,T]; Hq(0, 1)) (which is 
given by u(t, x) = g(x)u(t, x) in our case) can be treated with the same methods as 
the the homogeneous equation. It is also not difficult to verify that the map $ (and, 
by extension, is continuous and Frechet differentiable. We omit the technical 
details. Note that 

(4a) $(0^,0) = ^(T) = e^ tt(0„, e^ T M , 0) = 

(4b) 9 u *(^, e ^,0)-/i = 9 u $(0 M ,O)-ft (V/ieC([0,T];^(0,l))). 

Moreover, the derivative dJSf^, 0) : C([0, T]; f#(0, 1)) -> ^(0, 1) of $ at the point 
(0m, 0) is given by 

d,M^Q)-h = z(T;h) 1 
where 1 1— > z(t) = /i) is the solution of IB VP 

(5a) izt = -z xx + f(\Lp^\ 2 )z + 2/'(|^| 2 )Re(v? M ^)v9 + g(x)h(t,x) 

(5b) z(t,0) = z(t,l) = 

(5c) z(0,x) = 

3 This assumption can be made without loss of generality. To see this, assume that Theorem ^ 
is proved for this case. If g n , is the "actual" (non-smooth) indicator function for the interval f2', 
choose an open subinterval Q C Q' and a smooth function g with supp(g) C f2 and < g(x) < 1. 
Then, by Theorem ^ there will be a control u that solves the control problem (|la|l - ljldjl with £1' 
replaced by Q. Now u(t, x) := g(x)u(t, x) will be a suitable control for the original problem. 
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As usual, the time dependence of the term involving Re(...) is eliminated by the 
transformation 

z(t) := e~ itxt z{t), h(t, x) := e~ ilJt h(t, x), 
which, for f(s) = —s (cubic focusing nonlinearity) , gives the IBVP 

(6a) iz t = -z xx + fj,z- 2 2 - 20 2 Re(2) + g(x)h(t, x) 

(6b) z(t,0) = 2(t, 1)=0 

(6c) 2(0, a:) = 

We will work with the real (2 x 2)-system arising from (|6a |) -([Bc f) by decomposition 
in real and imaginary parts. (We will drop the // subscripts whenever there is no 
danger of ambiguity.) Consider the matrix operator 



(7) L 



-A + [i- cf) 2 (x) \_ ( L_ 

A-/i + 30 2 (x) J ~ \ -L+ 



this operator is to be understood as an operator in the Hilbert space H 1 x H 1 = 
[^(0,1)]* x [£#(0,1)]* with domain ^(0,1) x H%(0,1) = [H%(0, l)] 2 . Then eq. 
ml) takes the form 

lm(h(t, x)) 
-Re(h(t,x)) 



(8) Z t = LZ + g(x) 

where 



Z(t,x) 



'Re(z(t,x)Y 
Jm(2(t,x)) / 

The corresponding linearized control problem to be considered therefore reads. 



(9a) Z t = LZ + g(x)H(t,x) 

(9b) Z(t,0) = Z(t,l) = Q 

(9c) Z(0,x) = 

(9d) Z(T,x) = Z x {x) 

Our main objective is to prove the following controllability result for this linear 
system. 

Theorem 2 The control problem $9a \) - {Uc ty has a solution H = H(t,x), i.e. for 
any Z\ £ [i/^O,!)] 2 there exists a control function H £ C([0,T]; [Hq(0, l)] 2 ) such 
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that the solution Z G C([0,T]; [H%(0, l)] 2 ) of $a$)-$3) satisfies Moreover, the 

control H satisfies the estimate 

(10) sup WH^WhiKCWZ^. 

0<t<T 

Before proceeding with the proof of this theorem, we describe how Theorem^follows 
from it. 



2.2 Proof of Theorem Q] 

As mentioned in the introduction, our plan is to apply the IFT (see, e.g. [TBT Theo- 
rem 4.B]) to the map * : [H^O, l)] 2 x C([0, T}; H^O, 1)) -> #1(0, 1), *(y Q) y 1; u) := 
$(yo,«) - Vi- Since 

(11) e- iAiT 9 u $((/), 0) • e^ T h = z(T- 0, h). 

the map $„$((/>, 0) is onto iff the equation ([Sa j) -([Sc |) is exactly controllable, which, 
in turn, is the case precisely if the equation (J9aj) - (J9cjl is exactly controllable. Thus, 
by Th. refP, 5 U $(0, 0) (and hence d u ^((f), T ,0)) is onto. Moreover, 

119^(0,0) • h\\ = \\z(T)\\ m = \\Z(T)\\ H i > C\\H\\ m = C\\h\\ Hh 

which shows that the map d u <&(((>,0) is one-to-one as well and that its inverse is 
bounded. Thus, by the IFT, there exist neighbourhoods Y x Y 1 C [i^Q(0,l)] 2 and 
U C C([0, T\; H x (0, 1)) of (</>, e^ T </>) and 0, respectively, such that 

V(2/ , Vi) eY xY 1 3\u = u(y , Vl ) G U : V(y Q , y u u) = 0, 

which concludes the proof of the theorem (assuming that Th. |2]has been proven), 
since y\ G Hq(0, 1) is a reachable state iff there exists a function u G C([0, T]; H l {f), 1)) 
such that V(y , y u u) = $(y , u) - yi = 0. □ 



3 Exact controllability of linearized NLS 

The purpose of this section is to prove Theorem |21 The proof will employ the 
Hilbert Uniqueness Method (HUM) due to J.-L. Lions 0. We will be using the 
decomposition L = + ZA 1 ) of L, given by 
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We will also need the adjoint operator 

/ A-^ + 3(f) 2 (x) \ _ ( -L+ \ 

[iZ) L ~ V -A + M -0 2 (x) J~\L- J' 

The set-up is similar to the one in jB]. We use capital letters to remind ourselves 
that we are dealing with two-component functions. All equations below are to be 
understood subject to zero boundary conditions (i.e. are to be interpreted in the 
Hilbert space [Hq(0, l)] 2 ). The equation 

n - (W t = LW + g-V 

{ 6> \ W{T,x) = Z x {x) ' 

is decomposed into two "semi-homogeneous" equations: 



(14a) W t {1) = LW (1) , W {1) {T,x) = Z x (x) 

(14b) W t (2) = LW {2) +g-V } W {2) (T } x) = 0, where 

(14c) W = W (1) + W {2) 

The "HUM operator" is defined as 

(15) SV := -W {2 \0) 
where is the solution of (|14b|) with V given by 

(16) V t = -L*V, V(0) = V Q . 



If S possesses a (bounded) inverse (in the space [Hq(0, l)] 2 ), then the equation 

SV = W {1) (0) 

has a unique solution V , and W = W {1) + W {2) will satisfy (USD with W(0) = 0. So 
setting 

Z=W, H = V 

will solve the control problem (|9a )) -((9d )) . Moreover, the estimate ()10|) will follow 
from the apriori estimates (|27a|) and ()27b|) listed in Section 13.2. 1[ The crux of the 
HUM, therefore, consists in showing that S has a bounded inverse. This is done by 
showing that S satisfies an "observability estimate" of the form 

(SV , V ) > C\\V \\ 2 (VU G [HfcO, l)] 2 ) 

(w.r.t. the appropriate inner product and norm), which, by the Lax-Milgram 
Theorem, implies that S is an isomorphism of [Hq(0, l)] 2 . 
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3.1 The L 2 — observability estimate 

The space [L 2 (0, l)] 2 = L((0,1);C 2 ) is equipped with the standard inner product 
(., .), given by 

= ( U W , W (2) > L 2 ( (o,l),C) + V 2) >L2((0,1),C) 

l yl 

-u < - 1 - ) (x)'u( 2 )(a;)(ia; + / f ( - 1 - ) (x)i'( 2 )(a;)(ix; 

the corresponding norm is ||V|| = a/ (V, V). The objective is to show that there 
exists a constant Chum such that 

(17) (SV ,Vo) >C HUM \\V \\ 2 (W e [L 2 (0,1)] 2 ) 

The proof of (|17p will make use of the following properties of the operator L*. 

3.1.1 Spectral decomposition of — L* and solutions to (|16|) 

The justification of the properties listed in this section will be deferred to the ap- 
pendix (see IB.3jl . 

1. the spectrum of — L* consists of eigenvalues only 

2. all non-zero eigenvalues are purely imaginary 

3. all but finitely many non-zero eigenvalues are simple; there are no generalized 
eigenvectors associated with non-zero eigenvalues and each eigenspace is at 
most two-dimensional. 

Remark. In the sequel, we will for convenience assume that all non-zero eigen- 
values are simple 4 . All subsequent arguments (in particular, t^S] below) can 
easily be adapted to accommodate additional (linearly independent) eigenvec- 
tors which potentially occur for a finite number of eigenvalues. 

4. the multiplicity of the eigenvalue zero is 2; an eigenvector and a generalized 
eigenvector, V\ and W±, satisfying 

(18) -L*Vi = Q and - L*W X = V x , 

4 Thcre is persuasive numerical evidence to support this assumption. 
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are given by 

(19) v i= (t) and w i 



(in particular, V\,Wi G M 2 ). Moreover, these vectors form a basis of the 
generalized null space. 

5. Let V2, V3, . . . , V2, V3, . . . be the eigenvectors corresponding to the simple non- 
zero eigenvalues (see Remark in $03 above) A2, A3, . . ., A2, A3, Then Vi, W±, 

V2, V3, . . . , V2, V3, ... is a Schauder basis for [L 2 (0, l)] 2 as well as a Bessel se- 
quence, i.e. for any V G [L 2 (0, l)] 2 , there is a unique representation 

(20) V = C X V X + d X Wx + + ^V n ] ((Cn, C n ) G f) 

n>2 

and there exists a constant 5 > (independent of V") such that 

(21) | Cl |2+| rfl |2 + ^[| Cn |2 + 1^2] > £11^2 

n>2 

6. The (generalized) eigenf unctions Vi, Wi, V2, V3, . . . satisfy the uniform bound 



n>l 



22) m v := win < ini / g(x)\V n (x)\ 2 dx, / g(a:)|Wi(a;)pda: > > 



n ./n 



7. Let y G [L 2 (0, l)] 2 be real, V = c x V x + d x W x + E„> 2 [cnK + c n V n ] = Cl V x + 
d x W x + 2ReEn> 2 C «K and A « = i& (/?n £ K) • Then 



(23) n - 2 

= c 1 V 1 + td 1 W 1 + 2ReJ2c n e il3nt V n 



n>2 

is the unique solution of (fTE|). 
8. The sequence 

1, £, e iftli , e" iAlt (n > 2) 
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is a Riesz-Fischer sequence in L 2 (0,T; C), i.e., there exists a constant A > 
such that, for any £ 2 -sequence oi, b x , (a n , a n ) n > 2 C C, 



(24) 



/ a x + ht + l a n^ nt + ane-^ 1 } f dt 

J ° n>2 ' 

> A(\ ai \ 2 + |bi| 2 + ^[|a n | 



2 i I" |2' 



3.1.2 Proof of (HU) 

First we show that 

(25) 

Proof of H3J). Let 
Then 



T4EI.IT21 



S(t) := (^(t),^))- 
(^ (2) ,F) + (W« V t ) 

(LV^ 2 \F> + (gV,V) - (W®,L*V) = (gV,V) 



and so 



(SV ,V ) = -(W<V(0),V(0)) = (W^^,V(T)) - (W^(0),V(0)} 

=0 

= S(T) - S(0) = J* j t S{t)dt = j\gV, V)dt, 

which completes the proof of ()25|) . □ 

Now (fTTj) may be verified using standard arguments. The solution V(t) of ()27a|) has 
the representation 

V(t, x) = c x V x {x) + td x W x (x) + ^[c n e*V n (x) + ^e^V^x)] 

n>2 

(cf. (j2HJ))- We are going to apply with 

ai = CiVi(a;), b 1 = d 1 W 1 (x), a n = c n V n (x), d n = c n V n (x); 
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the result is 

(SV ,V ) ^ [ T (gV,V)dt= [ g{x) f \V(t, x)\ 2 dtdx 
Jo Jn Jo 

> A [ < 7 (x)(|c 1 | 2 |V r 1 (a;)| 2 + |d 1 | 2 |iy 1 (x)| 2 + 2^|c ri | 2 |K(x)| 2 )dx 

n>2 

> Am v i\c 1 \ + \d 1 \ 2 + 2 2_^\c n \ 2 \ > A B m v \\V \\ 2 , 

n>2 

which completes the proof of (|17|) . □ 
3.2 The ^-observability estimate 

Now that the L 2 -observability estimate (fT7j) has been established, it is sufficient to 
show that there exist constants C\ and C2 such that 

(26) ((SV ) x ,V 0x ) > CillVb^ll 2 - C 2 ||Vo|| 2 

Indeed this will imply 

((SV ) x ,V 0x ) > CMJ 2 - ^-(SV ,V ) 

_^ ^ 2 / ct r \r\ , nc^r\ t/ \ \ n in/ 112 



(SV , V ) + ((SV ) x ,V 0x ) ><7i||V 0ai | 



Chum 

The left-hand side is equivalent to the "natural" H 1 inner product 

(.,.) + <(•)*,(•)*>• 

Before presenting the proof of (|26j) we list the apriori estimates required. 
3.2.1 Apriori estimates 

We are going to need apriori estimates for the various functions involved. The proofs 
are standard fare and will be omitted. 

(27a) ||V(t)|| < C||Vb|| 

(27b) \\V x (t)\\ < C7(||V || + 11^0,11) 

(27c) ||^ (2) WII < C\\V \\ 

(27d) \\Wi 2 \t)\\ < C(||V || + 11^11) 
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(2) 

We will also need the equations for V x and W x , which are given by 



(28a) = LW&> + 4%< 2 ) + g x V + gV x , W®(T) = 

(28b) V xt = -L*V X -[LW]*V, V x (0) = V 0x ., 

3.2.2 Proof of (1251) 

We first prove an "if 1 analogy" of the identity 1)250: 

rT 



(29) 



((SVo)*,^) 



Proo/ o/ Let 
Then 



+ [ (g x V,V x )dt+ [ (gV x ,V x )dt 
Jo Jo 



S(t) := (WP(t),V x (t)). 



dt 



d ~ 



dt 



S(t) = (W^ V V X ) + (wp,v xt ) 



(LWP,V X ) + (L^W^,V X ) + (g x V,V x ) + (gV x , V x ) 
-(LWP,V X ) - (L^Wi 2 \V) 

(LUW®,V X ) - (L^Wi 2 \V) + (g x V,V x ) + (gV m ,V m ), 



which implies 
Using the abbreviations 



(30a) 
(30b) 
(30c) 

we estimate 
(31) 



h 
h 
h 



(LU W W,V X ) - (L^\V) 



(g x V,V x )dt 



dt 



(gV x ,V x )dt 



((SV Q ) x ,V 0x ) > h-\h\-\I 2 \. 



□ 
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Now, 



I'll < 



3||(0 2 ).||oo / (\\WW\\-\\V x \\ + \\WW\\-\\V\\)dt 



(32) 



< c\\v \\ ■ (\\Vo\\ + \\v 0x \\) 

Young 

< e\\V 0x \\ 2 + C e \\V \\ 2 . 



The estimate for I 2 is even simpler. 



(33) 



r T I27al . l27bl 

\h\ < H&lloc / ||v|| • IIKIIrft < c\\v \\ ■ (11^.11 + livw! 

Jo 

Young 

< e||^,|| 2 + C £ ||^ || 2 



Finally, we need to estimate J3 from below. To do this, we write the solution of 
(I28bl) in Duhamel form 



(34) 



V x (t) = e~ LH V 0x - 1 < 



-L*(t-s) 



[LU]*V(s)ds, 



where e L *' denotes the semi-group corresponding to the equation If we set 

V (0 \t) := e~ LH Vo x , 

we get from (fT7|) and (|23j) that 

(35) [ T [ g(x)\V^(t,x)\ 2 dxdt = [ T (gV^\V^)dt > C HUM \\V 0x \\ 2 



J n 



(here we mean by |.| the norm of C 2 , i.e. I (") I = \u\ 2 + |f| 2 .) Now, 



|/ 3 | > f f g(x)\V w (t,x)\ 2 dxdt- 
Jo Jn 




03J 

> CWmIIVqzI 



'0 Jn 

r-T 



-L*(t-s)rr(l)]* 




Jn 



[L^)*V(s,x)ds 



-(ts)[ L U}*V(s,x)di 



dxdt 



dxdt 
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To finish the proof, we need to estimate the integral term. 

-T r rT pT 



[ [ g{x)\...\ 2 dxdt < [ [ [ |e- L *(*- s )[LW]*y(s,x)| 2 dxdsdt 

Jo Jn Jo Jo Jn 

= [ [ \\e- L * {t - s) [LW)*V(s,.)\\ 2 dsdt 

Jo Jo 



JO 

^ ° r / T ii[4 i) ]*^( s >-)ir^ 

Jo Jo 

rj-i ^^^^^^ 

< C\\(<P 2 ) X \\ 2 [ \\V(s)\\ 2 ds < c\\v \\ 2 
Jo 

Thus, 

(36) / 3 > C^ Af ||Vb x || 2 -C||Vo|| 2 

and so 

((SV ) x ,V 0x ) > h-\h\-\I 2 \ > (C HUM -2e)\\V 0x \\ 2 -C £ \\V \\ 2 , 
which will conclude the proof of (ffilj) if e is chosen sufficiently small. □ 

4 Concluding remarks 

There are a number of modifications/generalizations of the control problem (fTaj) - (jld|) 
which are of interest. Those include 

• space dimension > 1; 

• other boundary conditions such as periodic boundary conditions; 

• the ground state M in the assumption of Theorem ^ may be replaced with 
some excited state fsee IA.1|) ; 

• zero (or "box" ) boundary conditions may be interpreted as an infinite potential 
well, which one may want to replace with other, typically confining, potentials, 
such as the harmonic-oscillator potential (in this case, one would reasonably 
set fl = (— oo, oo).) 
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Any of these modifications will obviously change some of the spectral properties in 
13.1.11 To see whether a controllability result as in Theorem ^ can be proved under 
these modified assumptions as well, the ramifications for the application of HUM 
will have to be carefully examined. 

Appendix 

A Bound states 

We define bound states as real solutions of the BVP (j2a j) .(|2b jl . 

A.l Bound states in terms of elliptic functions 

It is well-known that explicit formulas for the solutions of ()2aj) . (j2hj) are available 
in terms of Jacobian elliptic functions; see, e.g. |2j. If j G {0, 1,2, . . .}, then (f>j(x) 
will denote the (real- valued) solution of (J2U , (J2EJ) which possesses precisely j zeros 
("nodes") within the interval (0,1). The node-less solution <fi := O is referred to 
as the ground state] the solutions <f>j {j > 1) with one or multiple nodes are called 
excited states. To find an explicit solution formula for <f>j, we first solve the equation 

(37) ^ = A{j + lf{2k 2 -l)K{kf 

for k, where K(k) denotes the complete elliptic integral of the first kind (see, e.g. 
PP). Note that, since K(k) is a strictly increasing continuous function of k G [0, 1) 
satisfying lim^i- K{k) = oo, equation (|3T|) has exactly one solution k = kj{fjL) 
for any choice of parameters /i > and j G {0, 1,2, . . .}. Moreover, the function 
kj : [0, oo) — > [-^=, 1) is continuous and strictly increasing as well, and satisfies 
lim^oo kj(s) = 1. Now the solution (f)j of (J2~aj) . (}2b"|) is given by {k = kj{jj)) 

,38) *w - ;fe m (ffcf 

= 2V2(j + l)kK(k) cn(2(j + l)K(k)(x - \) + \j] 2 K(k),k) , 
where [j}2 := j mod 2. The following properties are readily proved. 



15 



Lemma 3 Let ip £ Hq(0, 1) 6e a non-trivial weak solution of \2a\) . Mb]), i.e. 

[ [il> x v x + fj,il>v-if! 3 v]dx = (Vve£#(0,l)). 
Jo 

T/ien 

(%) -0 e C^O, 1) n C[0, 1], i.e. ?/> zs a classical solution of pta)). (flfcj) . 
(wj There exists j G {0, 1,2,.. .} sttc/i £/iai ^ /ias exactly j zeros in (0, 1). 
(raj ^(x) = <j>j{x), where <pj is given by tffify . 
We also need the following convexity property. 
Lemma 4 TTie function \x ^ \\4>\\2 zs an increasing function of \x. 
Proof. This follows from the identity 

/K(k) 
cn 2 {y,k)dy 
■K(k) 

and the fact that k = k(/i) is a strictly increasing function of /i. □ 

A. 2 Variational description of the ground state 

There are various variational descriptions of the ground state in the whole-space 
case; see e.g. O 4.2] and [3 8.1]; in the zero-boundary case we find the formulation 
presented in Lemma El below to be a useful one. Let fj,\ be the smallest eigenvalue 
of the ID Laplacian — on (0, 1) with zero boundary conditions, i.e. /ii = n 2 . 

Lemma 5 Let fi > —fit. Then there exists a positive minimizer 

w E{ue #o(0, 1) | f u 4 dx = 1} =: M 
Jo 

for the constrained minimization problem 

1 f 1 

inf - / [u x 2 + nu 2 ]dx, 
ueM 2 J 

and the unique positive solution (j) of |B?q }) . $2b lj) is given by 0(x) = A 1 / 2 w(x) ; for some 
suitable A > (Lagrange multiplier). 

Proof. See e.g. [11, Theorem 2.1 and proof]. □ 
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B Spectral properties of — L* 



B.l Properties of L + and L 



The operators 



(40a) 



dx 2 



+ /x- c/) 2 (x) u, w(0) = = 



(40b) 



L + u 



dx 2 



+ H — 3<fi 2 (x) u, u(0) 



u(l) = 



are regular Sturm-Liouville (SL) operators, so we can make use of the SL theory. 
Lemma 6 (i) ker(LJ) = span{(p} 
(ii) L_\[ span ^± > 

Proof, (i) Clearly, L_0 = 0, so is eigenfunction for L_ with eigenvalue A = 0. 
Since all eigenspaces are one-dimensional, the assertion follows. 

(ii). Since <p has no zeros in (0, 1), it is an eigenfunction for the smallest eigenvalue; 
hence Ai = and A n > (Vn > 2), which implies assertion (ii). □ 

Remark. If <fi is not the ground state but an excited state with j G N \ {0} nodes, 
then the operator L_ will have exactly j negative eigenvalues. 

Lemma 7 The operator L + has exactly one negative eigenvalue and all the other 
eigenvalues are positive. In particular, ker{L + ) = {0}. 

Proof. The proof proceeds in four steps. 

Step 1. L + possess at least one negative eigenvalue. This follows from 



and the minimax principle. 

Step 2. (L + r],r]) > for all rj e [span(0 3 )] ± . This is a slight adaptation of the 
arguments in ^QJ Sectionl3]; we use similar notation. Let the functionals J[u] and 
W[u] be defined by 



and let w G M = {u G Hq(0,1) \ W[u] = 1} be a positive minimizer of the 
constrained minimization problem inf^g^ J[u]. By Lemma El the ground state <p is 



(L +( f>,<f>) = (L_0,0) - 2\\m 



2|M|J<0 
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given by <j){x) = X 1 ^ 2 w{x) for some positive constant A, which arises as a Lagrange 
multiplier. Now let rj G [span(0 3 )] ± ; we write rj in the form rj = w := -j^\ z= ow(., z) 
where z i— > w(.,z) is a smooth curve in Z/q(0, 1) such that = 1 (i.e. 

w(., z) G M) for all z, and w(., 0) = w. Since w is a minimizer, we have 



(41) 
and 



dz 



2=0 



J[w(. : z)\ = I [w x w x + fiwwjdx 



< 



dz 2 



(42) 



2=0 
1 



J[w(.,z)) 



[w x + w x w x + fi(w + ww)]dx 



[(w xx - fiw)w + {w xx - /j,w)w]dx 



where w := -i^\z=ow(., z ) Moreover, from the constraint z)] = 1 we get 



(43) 
and 
(44) 



d 

dz 



2=0 



W[w(., z)] = / w wdx 







dz 2 



2=0 



W[w(.,z)} = I [3w 2 w 2 + w 3 w]dx. 



Also, the Lagrange-multiplier rule implies 

(45) 



/ [w x v x + fiwv] dx = X w 3 vdx (Vi> G Hq(0, 1)) 
Jo Jo 



-w xx + imw\w dx 



f .„\ v=w 

(46) 

Inserting this into ()42|) gives 
< 



X I w wdx 



X I 3w 2 w 2 dx 



-w xx + fiw — X3w 2 w]wdx 

-Vxx + m - 30 2 ri]ridx = (L + r/, rj) 



(Note that choosing v = w in (}45|) also yields J Q [w 2 + fiw 2 ] dx = X C w 4 dx implying 
that A > 0, since /i > —fii.) 
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Step 3. The second eigenvalue, A2, is non-negative. This can be shown by repeating 
word-by-word the proof in (TU1 page 58] if \l/ is defined as 3 and II is interpreted as 
the orthogonal projection onto the subspace [span(0 3 )] ± . 

Step 4. A2 > 0. Assume that A2 = and let v denote an eigenfunction for 
A2. By the symmetry of <p(x), we may assume w.l.o.g. that v is either odd (i.e. 
v(l — x) = —v(x)) or even (i.e. v(l — x) = v(x)). In the first case, we have = 
and v coincides with a constant multiple of 4>' by ODE uniqueness. However, this is 
impossible, since does not satisfy zero boundary conditions. We are therefore left 
with the second case (v even). Since v is an eigenfunction for the second eigenvalue, 
it has precisely one zero in (0, 1) by standard SL theory. By symmetry this zero 
must occur at x — ~, which is impossible as we saw above. □ 



B.2 Eigenvalues and eigenf unctions for n 



00 



The properties El and |H1 in Section f3. 1. II are based on asymptotic (n — > 00) formulas 
for the eigenvalues A n ,A n and eigenfunctions V n (x), V n (x) of — L*, which are given 
in Lemmas |H1 and El below. We are going to make use of the fact (proved in IB. 31 
below) that all non-zero eigenvalues are purely imaginary, i.e. we write 



(47) 



An ij3rif A r 



-if3 n , (3 n e R, p n > 0. 



Moreover, it will be convenient to employ a similarity transformation ^0J (12.15)]: 
Let 

' 1 1 



J := 



i —1 



Then 

(48) 
and so 



-iL* = J ■ 



A - fj, 






-A + /J, 



+ 



J- 1 . 



spec(-L*) = zspec(M), V n = JW+, V n = JW~ n 



where (±/3„, W^) are the eigenpairs for the operator 



M :-- 



A - fi 






-A + /i 



+ 



=: M(°'+I (1) . 



i.e. 



(49) 



MW± = ±P n W 1 



± 

n ' 
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Writing j3 = ±/?„, W = = ("), the characteristic equation is equivalent to 
the BVP 



(50a) 
(50b) 

Note that 
(51) 



u"-(n + (3)u = -(p 2 (2u-v), u(0)=u(l) = 
v"-([x-P)v = <p 2 (u-2v), v(0)=v(l) = 0. 



w: 



w: 



Lemma 8 (i) The operator M (and hence —L*) is a spectral operator. More 
precisely, the collection of eigenvectors and generalized eigenvectors for M 
forms a Schauder basis for L 2 (0, 1;C) and all eigenvalues with sufficiently 
large indices n are simple. 

(ii) There exists an index n G N and a constant C > such that 

KnV + ix) - /3 n \ < C 

for all n > n . 

Proof. The operator is self-adjoint; its eigenvalues and eigenfunctions are given 
by 



sm(nnx)(° i ), W^ Y 



sm mnr 



(where M^w!i° ]± = fi^W^). Thus, the operator M = M< ' + M« "is" a 
bounded perturbation of a self-adjoint operator whose spectrum consists of simple 
eigenvalues only. The assertion now follows from [HJ Th. 4. 15. a]; see also [S]. □ 

Lemma 9 There exist an index no G N and a constant C > such that f3 n > n and 



(52a) 
(52b) 



W+(x)-sin (Va^z) (J) 
W~[x) - sin (y@ n - ix xj fjj 



< 



< 



c 

n 
C 
n 



for all n > n . 
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Remark. The assertion of the lemma may be expressed in a more intuitive, if 
slightly informal, manner by means of the asymptotic formulas 

W+{x) = sin(u; n x) (J) + 0{n- 1 ), n -> oo 
W~(x) = sin (oj n x) Q + 0(n _1 ), n — > oo. 

(uniformly in x and n) where u; n := y/ (3 n — ~ 7m, as n — > oo, by Lemma |H1 (ii). 

Proof. Clearly, by Lemma[S](ii), lim n ^oo j3 n = oo; so we may assume that (3 n — jx > 0. 
Let := y/ f3 n ± fx. Because of (JoTjl it is sufficient to consider W+; write = 
("r ). Viewing the R.H.S.'s of (|50aJ) . (j50b|) as inhomogeneities, we write the system 
as 

(53a) u" - [u)+fu = f(x), u(0) = u(l) = 

(53b) v" + [uj-} 2 v = g(x), v(0) = v(l) = 0. 

where 

(53c) f(x) = (p 2 (x)[v(x) — 2u(x)} and g(x) = <p 2 (x){u(x) —2v(x)]. 

Note that the homogeneous BVP associated with eq. ()53b|) has the solution sin(o;~a;), 
while the homogeneous BVP associated with eq. ()53a)l does not have any non-trivial 
solution, which implies that there is a Green's function, T u) +(x,^), associated with 
eq. ()53aj) . Utilizing this Green's function and the Duhamel Principle, solutions to 
(|53a|h(|53b[) may be written as 

(54a) u(x) = I r^(x,0mH 
Jo 

i r 

(54b) v(x) = csm(uj n x)-\ / sin(u; n (x - OMO 6 ^ eel) 

u n Jo 

By the linearity of the system f|5Ua|) . f)50b|l . we may assume that c = 1. Thus, 

(55a) u(x) = I r ui (x,0A0[M0-^M 

Jo 

(55b) v(x) = S in(u-x) + — [ sm(Lj~(x-0)(f) 2 (0H0-M0n 

w n Jo 
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Standard calculations yield the explicit formula for the Green's function T w + , which 
is given by 



(56) < 



& = a + \( + , {[sinh(a;+Oe-^ + - sinh(u,+(l - £))]e<^ 

-[smh(aj+£)e w " - sinh(o;+(l - 0)]e" w " 
sinh(u;+|a; — £|) 



2u: 



The function |r w +(x,£)| assumes its maximum on [0, l] 2 at the point (x,£) = (|, |) 
and its maximum value is given by 



(57) lr + f 1 Ml = sinh2( ^ } = cosh ^" + ) ~ 1 =o(—\ ( U - 

1 J l 1 "*^'^! ^+sinh(w+) 2u;+ sinh(u;+) \ -+ ^ lr ' 



a;. 



oo 



(uniformly in o;+). Now it is a matter of routine estimates (combining the repre- 
sentation ([55a|h(|55bp with the uniform estimate (jBTj) and the property uj^ ~ n, as 
n — > oo) to verify the assertion of the lemma. □ 

B.3 Verification of the properties listed in Section 13.1.11 

n Lemma |Hl(i)- 

[21 Let A £ C \ {0} be an eigenvalue for L. Then the argument on page 13 of 
[3] ("Claim") shows that A 2 is real. We can then use the argument on pages 
49 and 50 of ^U] to show that A itself is purely imaginary. This argument 
uses, besides the properties of L + listed in Lemma the convexity property 



_d_ 



2(<9 M 0,0) > (Lemma El). 



03 The first assertion (all but finitely many eigenvalues are simple) was already 
mentioned in Lemma |Hl(i). 

A proof of the second one (there are no generalized eigenvectors associated 
with non-zero eigenvalues) can be found in [TO! pages 50-51]. The proof refers 
to the whole-space case, but carries over to the zero-boundary case if the 
domains of the various operators involved are modified suitably. 
The third assertion (each eigenspace is at most two-dimensional) may be seen 
as follows. Let A £ C be an eigenvalue and Eig(A) the corresponding eigenspace 
and consider the linear map F : Eig(A) — ► C 2 , defined by FV := V'(0). Using 
the left boundary condition (V(0) = Q) and ODE uniqueness, is is easy to 
see that F is one-to-one. Thus dim(Eig(A)) = dim(im(F)) < 2. 
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EJ Let (") be an eigenfunction for L* associated with the eigenvalue A = 0. Then 
L_u = and L + v = 0. From Lemma|ni(i) and Lemma[7|we get u = ccj) (c G C) 
and v — 0; hence 

") e span { (o 

It is easy to see that W\ = satisfies — L*W\ = V\\ it is therefore a 

generalized eigenvector. We want to show that there cannot be another linearly 
independent generalized eigenvector. To prove this, let's assume that (") is 
such an "additional" generalized eigenvector, i.e. 

L W Wl 

since the eigenspace is one- dimensional. In particular, 

(58) L_u = -d^. 
By the Fredholm alternative, this implies 

(59) O = 2(0,<V/>> = !-(0,0> 

since is a non-trivial solution to L_</> = 0. But this contradicts Lemma 0] As 
a result, no solution u to (|5Sj) can exist, so there is no "additional" generalized 
eigenvector. 

Basis property: Lemma |Hl(i). 

The Bessel-sequence property follows from Lemma 03 Since it is obviously 
sufficient to establish that the sequence of (generalized) eigenvectors for M is 
a Bessel sequence, we will show that 



(60) W = □ G [L 2 (0, l)] 2 : (ai) nm := ((W, W±» n > N G 



Clearly, we may skip a finite number of terms in (|60|). For simplicity, we will 
also restrict ourselves to a+. Let n be an index such that 

W+{x) = sinKx) (°) + (Vx G [0, 1], n > n ), 
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where u n := \f(3 n — ju and i?+ : [0, 1] — > C 2 are continuous functions satisfying 
|-Rn( x )l — C uniformly in x and n (see Lemma Now 

a n = v \ x ) sm{u; n x)ax H =: o„ + c„ . 

./,) n 

The sequence (c+) is easily seen to be an ^-sequence: 

\(W,R+)\<C I \W{x)\dx<C\\W\\ \4\<^^e£ 2 . 



n 



To see that (6+) is square-summable as well, we note that, by Lemma |H1 the 
sequence (uj n ) n > m has the asymptotics u) n ~ 7m, as n — > oo, and is therefore 
separated. By jTHl Theorem 3.4], this implies that the exponential system 
{e luJnX } n > no forms a Bessel sequence in L 2 (0, 1). Thus 

b+ n := I Vj (x)e^ dxef (JG{1,2}), 
Jo 

where V\{x) := Re(t>(x)) and v 2 (x) := lm(v(x)), and so 

&+ = Im(&+J+*Im(&+ n )G£ 2 , 

which completes the proof of (a+) G £ 2 . 

IH1 This also follows from Lemma El Let Q' = (a, b) C (0, 1). Then we have 

f ■ 2/ \ j & - « 

/ sin [uixjdx — > , asw-> oo. 

Jn> 2 
As above, we appeal to Lemma El to write in the form 

W+(x) = sin(u; n x) Q + (Vx G [0, 1], n > n ). 

Now 

|W+(x)| 2 dx > 



sm 2 (u n x)dx / \sm(u} n x) ((°!) , R^(x))\dx — 

nj n > 

-- 2 [ \R±(x)\ 2 dx 

f 2 , x , 2C(b-a) C 2 (b-a) b-a 
> / sin 2 (u; n x)rfx - ■ 



n n 2 2 

as n — > oo (note that lim^ooU^ = oo). This concludes the proof. 
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clear 

IH1 According to Lemma |HJ we may write 

P n = 7r 2 n 2 + r n (n>n ), 
where \r n \ < C for some constant C > (independent of n). It follows that 

fin+1 ~ fin = 7T 2 [0 + l) 2 - U 2 } + T n+1 - T n = 7T 2 [2n + 1] + r n+1 - T n -> OO, 

as n — > oo, which, by jU Corollary], implies that the sequence 



g ■= { e - iftl *,l,e^} 



n>2 



is a Riesz-Fischer sequence in L 2 (0,T) for every T > 0. We will show that 
adding the function f(t) = t to S will result in a Riesz-Fischer sequence in 
L 2 (0,T) as well, ifT>0 is large enough. For ease of notation, let's define 

e (t) = l, e ±m (t) := e ±i/3 ™ +1 * (m > 1). 

The fact that the sequence <f is a Riesz-Fischer means that there is a constant 
A/g > such that 



(61) 
Thus 



T 



^ a n e n (t) + 6/(t) 



dt> Ag |a„| 



ED 



(V(a. 



dt > ^||a|| 2 2 + ||/||i 2(0iT) |6| i -2|&|||a|| £ 2||F|| £2 



where the sequence (-F n ) ngZ is defined by 

-T />T 



e±m(t)f(t)dt = I ( 

</0 



±i0m + lt 



f(t)dt (m > 0). 



(Since the sequence (— /? m , /3 m ) me N is obviously separated, the sequence ^ is 
a Bessel sequence, which implies that F G £ 2 . Since f(t) = t, the F n 's can 
also be computed explicitly to verify the square-summability of F.) Utilizing 
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the elementary Young's inequality, we can continue the estimation above to 
obtain 



T 



J2 a n e n {t) + bf{t) 



dt > e||a||£2 + 



L 2 (0,T) 



A,, 



\b\ s 



(for e G (0, A?)) j which will yield the assertion provided that the condition 



(62) 



|2 

I P 



<{A g -e 



2 

L 2 (0,T)- 



is satisfied. Clearly, 



since f(t) = t. Moreover, it is a matter 



12 

lL 2 (0,T) 3 ' 

of routine calculations to verify that < CfT 2 for some T-independent 

constant Cf. Thus, condition ()62|) takes the form 



(63) 



C f <%(A,-e) 



(note that the constant Ag depends on T, Ag = Ag{T)). Finally, it is easy 
to see that Ag can be chosen such that AgiT^) > Ag(T\) if Ti > 2\, which 
implies that condition (jfjHjl (and hence (|62|l) can be fulfilled by choosing T 
large enough. This concludes the proof that the sequence $ U {/} is a the 
Riesz-Fischer sequence in L 2 (0,T) if T > is sufficiently large. 
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